Vortex core transitions in superfluid 3 He in globally anisotropic aerogels 
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Core structures of a single vortex in A-like and B-like phases of superfluid 3 He in uniaxially 
compressed and stretched aerogels are studied by numerically solving Ginzburg-Landau equations 
derived microscopically. It is found that, although any uniaxial deformation leads to a wider A-like 
phase with the axial pairing in the pressure-temperature phase diagram, the vortex core states in 
the two phases in aerogel depend highly on the type of deformation. In a compressed aerogel, the 
first-order vortex core transition (VCT) previously seen in the bulk B phase appears at any pressure 
in the B-like phase while no strange vortex core is expected in the corresponding A-like phase. By 
contrast, in a stretched aerogel, the VCT in the B-like phase is lost while another VCT is expected 
to occur between a nonunitary core and a polar one in the A-like phase. Experimental search for 
these results is hoped to understand correlation between superfluid 3 He and aerogel structure. 



I. INTRODUCTION 

Superfluid 3 He is in a spin triplet Cooper pairing state 
with a p-wave orbital symmetry 1]. Recent research 
interests have been directed toward superfluid 3 He in 
an aerogel with highly porous structure composed of 
silica strands to search for possible impurity-induced 
or anisotropy-induced novel superfluid states 043- At 
least, two superfluid states, A-like and B-like phases, 
seem to appear in aerogel and their pairing states are 
the same as those of the bulk A and B phases, i.e., the 
axial 0-0] and isotropic pairings, respectively. On the 
other hand, by introducing some anisotropy over large 
scales in aerogel, a change in pairing state might occur. 

In aerogel, locally anisotropic scattering events due to 
silica strands may stabilize the A-like phase, although 
this effect is mostly cancelled by an impurity-induced re- 
duction in the strong-coupling (SC) correction |g, Is 
necessary for stabilizing the axial pairing. Recent stud- 
ies have demonstrated that a global anisotropy existing 
over large scales helps to stabilize the A-like phase be- 
cause it promotes an orientational ordering of the orbital 
axis 1 of the axial pairing state and may screen quenched 
disorder effects of aerogel on 1. For instance, in a uniax- 
ially compressed aerogel, the 1 vector of the A-like phase 
tends to be aligned along the compressed direction while, 
in a uniaxially stretched aerogel, 1 should lie in the plane 
perpendicular to the stretched direction. Such a control 
of the 1 texture due to a deformation seems to have been 
verified through recent NMR measurements [ljj, [ll| . In 
addition, as the stretched deformation is enhanced, the 
novel polar pairing state may be realized in a narrow 
temperature region just below the superfluid transition 
temperature T C {P) [5|. Although the polar pairing state 
has not been detected yet experimentally, a wider A-like 
phase and an emergence of the polar pairing state might 
be indicated in spatially local regions as a consequence 
of the multiplicity of components of the order parameter 
in 3 He. 

In this work, core states of an isolated vortex in super- 
fluid 3 He in globally anisotropic aerogels are examined 



based on a Ginzburg-Landau (GL) free-energy functional 
Hgl derived microscopically. By assuming a boundary 
condition consistent with the radially symmetric core- 
less vortex [12| , we find that a polar-core vortex and 
the resulting vortex core transition (VCT) in the A-like 
phase are obtained in a uniaxially stretched aerogel. In 
contrast, in the A-like phase in uniaxially compressed 
aerogel, the 1 vector tends to align along the anisotropy 
axis so that no intriguing vortex core structure results in. 
Therefore, regarding core structures in the A-like phase, 
we will limit ourselves to possible situations in the A-like 
phase in uniaxially stretched aerogel. We also examine 
possible vortex core structures in the B-like phase and 
find that the VCT found in the bulk B phase [H, [llj| . 
which occurs between a nonaxisymmctric structure and 
an axisymmetric one with A-phase core, survives in a uni- 
axially compressed aerogel even in the low-pressure limit 
while such a VCT accompanied by the A-phase core is 
not realized in the stretched aerogel. 



II. GINZBURG-LANDAU THEORY 

Our starting point is the weak-coupling BCS Hamilto- 
nian with a p-wave pairing interaction term and a term 
describing impurity scatterings brought by the aerogel 
structure [3|, |5[ . The latter term is expressed in the sim- 
ple nonmagnetic and potential scattering form 



ffimp = / d 3 ru(r)h(r), 



(1) 



where n is the quasiparticle number density. The scat- 
tering potential u has zero mean and yields the ensemble 



\u k \i = u z (l + 5 u (k ■ zY) 



(2) 



at a fixed k, where u^ is the Fourier transform of u(r), 
k = h/kp, the overbar denotes the random average, 
and a small parameter S u measuring the nature of a 
global anisotropy is negative (positive) for uniaxially 



stretched (compressed) aerogels [5J. Within the sim- 
plest Born approximation, the relaxation rate becomes 
T~ l = 2ttN(Q)uI(1 + 5 u /3), where N(Q) is the density 
of state per spin on the Fermi surface. Although, as 
a model for real aerogels, Mk should have more compli- 
cacted k dependence, other k dependence expressing the 
local anisotropy of aerogel will be neglected hereafter by 
assuming the global anisotropy to screen the quenched 
disorder effects of aerogel on 1. The GL free-energy func- 
tional "Hql has been derived elsewhere [5[ from the BCS 
Hamiltonian with the anisotropic impurity scattering po- 
tential defined by Eq. @ and is expressed as 
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, (3) 



where A^j is the order-parameter field with a spin in- 
dex 77 and an orbital one i which denote the Cartesian 
coordinate axes, 
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where T c q is the superfluid transition temperature in 
bulk, d T = 1/(4tiTt), /9b (T) = 7C(3)iV(0)/(2407r 2 T 2 ), w F 
is the Fermi velocity, ip^ n '(z) is the nth order digamma 
function, and ip(z) = tp^°'(z). In Eq. @, any term de- 
scribing quenched disorder effect on A n j of the locally 
anisotropic aerogel has been dropped by favoring effects 
of a global anisotropy of aerogel. Further, higher order 
contributions in dr 5 U are also dropped in Hgl by as- 
suming them to be negligibly small. 

The SC correction 5f3j to /3j needs to be incorporated. 
Effects of an aerogel on 5(3 j have been studied thoroughly 
in Ref . \G\ by examining a general form of Sf3j , and it is 



found that the impurity effect other than the relaxation 
rate t _1 , which is a contribution to the quasiparticle ef- 
fective interaction ja], leads to a reduction in \6f3j\. Here- 
after, the expression of S/3j in the spin-fluctuation model 
15] will be used for simplicity. In this model, we have 
Sfc = (6p 2 + 56Pi)/6, <S/3 4 = Sfo + 5601, and 5f3 5 = 78fa, 
and 

6fa = -3.3 x 10- 3 5/3 (T)^^[(M d) (m)] 2 , 
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where (5 is defined in Eq. (3.10) of Ref. [l5| and is scaled 
with T c o/Ef and Ep is the Fermi energy. Then, we ex- 
press S = r]T c o/Ep and take n as a parameter measuring 
the strength of the SC correction. Below, the value 300 
will be used as 77 of the bulk liquid together with T c o(P) 
(Ref. [1]) and E F (P) (Ref. 0) data. In the aerogel 
case, a smaller rj value will be assumed to mimic the re- 
duction in the SC correction in aerogel Q. 

Below, we study the structure of an isolated vortex par- 
allel to the z axis by solving the GL differential equations 
under a suitable boundary condition. Following Salomaa 
and Volovik's approach [14| , we solve the equations in the 
case with an axisymmetric vortex and then, determine 
possible vortex structures by introducing a nonaxisym- 
metric part of each A^^ as a perturbation and examining 
the stability of the axisymmetric solution against nonax- 
isymmetric ones. To proceed, it is convenient to rewrite 
%gl by representing A v ^(y) in terms of the cylindrical 
frame as follows: 

A Vti (r) = A{T)J2 J2 ^ A i^ 1 ( r ) eig ' ei(m """*. 
Q W 

(6) 
where A(T) is the energy-gap amplitude far from the 
vortex axis, m is the circulation quantum number, Q 
is an integer, /1 and v specify the projections of the 
spin and the orbital angular momentum of the Cooper 
pair, respectively, and take three values (+1, 0, -1) with 
A? = [xi ± iyi)/s/2, A" = Z{. Note that, hereafter, all 

components of CJ^fJ (r) have only radial dependences and 
that axisymmetric and nonaxisymmctric parts of A rh i 

correspond to CJifJ(r) with Q — and with nonzero 



Q, respectively. Equation ([6]) shows that the axisym- 
metric component C^l with m — /i — v 7^ does not 
become nonzero at the vortex center where r = since <p 
becomes multi- valued there while the components with 
m — [i — v = can exist at r = 0. In contrast to the case 
of s-wave superfluid where a vortex always has a normal 
core, an axisymmetric vortex in superfluid 3 He can have 
nonvanishing superfluid components with m — /i — v = 
at the vortex center instead of the normal core so as to 
gain the condensation energy. Hereafter, a vortex whose 
center is occupied with a superfluid component X dif- 
fering from the components realized far from the vortex 
center will be called as "a X-core vortex." 

In rewriting %gl in terms of Eq. ((6J, nonaxisymmctric 
(Q ^ 0) components C^fJ with the same \Q\ couple with 
each other, while those with different \Q\ are not mixed 



at the lowest order in the correction C, 



(<2#0) 
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rewrite the GL free-energy functional per unit length in 



the 2-direction Hq L as follows: 
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|Q| < 2 + m so that higher \Q\ components are ne- 
glected. Below, we show results of our calculation for 
m = 1 vortices in deformed B-like phases and the uniax- 
ially stretched A-like phase by assuming each of them to 
be an vortex in which all components of Cjffj are real. 
In obtaining the results, the parameters 1/(2ttt) = 0.135 
(mK) and rj = 266 are used in the aerogel case, and the 
vortex radius f& — 20 is assumed. 



III. VORTEX CORE TRANSITION IN B-LIKE 
PHASES IN ANISOTROPIC AEROGEL 



First, a possible VCT in the B-like phase in uniaxially 
deformed aerogels will be considered. The uniaxial de- 
formation also affects A Vt i far from the vortex core, and 
thus, the boundary condition we choose on A n ^ far from 
the vortex core is 
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keeping the normalization 20 2 

mined by minimizing the free energy in the vortex-free 

situation. 

In the bulk B phase where 1/r = and O xy = O z = 
l/v3, we obtain two solutions of a single vortex, an 
axisymmetric vortex and a nonaxisymmetric one with 
\Q\ = 2 components, in agreement with the previous re- 
sults [I3, [l4[ • The former vortex gains a large | Cg _j_ | and a 

small I C + Q I at the vortex center which correspond to the 
A-phase component with 1 || z and the /3-phase 1| one, re- 
spectively. Thus, this structure is called an A-phase-core 
vortex. On the other hand, the latter vortex has nonva- 



nishing four components Cq_j_, Cq_ , 



and (M 



r (-2) 



where summations for g, are carried out under the con- 
straint qi + q 2 - q-3 - <?4 = 0, a„ = a n /\a + a z 2 z \ 
(n = 0,z), A = |A(T)| 2 A/|a +a z O% f = r/£ GL (T), 
^gl(T) = {K/\a + O 2 z a z \) 1 / 2 , and \O z \ is the magnitude 
of A ViZ far from the vortex center (see below). 

We set the boundary condition at fj, far from the vor- 
tex center and numerically solve the differential equations 
obtained by varying Eq.fJT]) with respect to CJ^fJ . As in- 
dicated in Eq. ([7]), we can examine the instability of an 
axisymmetric (Q — 0) solution against nonaxisymmet- 
ric corrections Cjffv ' at each \Q\ value independently 
Furthermore, any component with \Q\ > 3 + m cannot 
survive at the vortex center since <f> becomes multi- valued 
(±|Q| + m — /i — v 7^ 0) there. Since we are interested 
in the order-parameter structure close to the vortex cen- 
ter, for brevity, we will incorporate only corrections with 



«at the vortex center, forming a twofold symmetric vor- 
)tex core state, and thus, is called a double-core vortex. 
''We also compare the total free energy J dff /(f) of the 
Hxisymmetric vortex with that of the nonaxisymmetric 
vortex and find that the axisymmetric A-phase-core vor- 
tex is stable at higher temperatures. The resulting tran- 
sition (VCT) curve between the two vortex core states 
is denoted by a thin (black) solid curve in Fig. 1. The 
stability of the A-phase-core vortex against the double- 
core vortex results from the large gain in the conden- 
sation energy of its core phase overcoming the energy 
cost due to the spacial variation in the order parame- 
ter larger than that for the double-core vortex [l3j. The 
solid black curve in Fig. 1 remains nearly straight, re- 
flecting the bulk A phase stabilized by the SC correction, 
namely, the pressure dependence of T c q/ ' Ep in 5 in the 
spin-fluctuation model. However, the experimentally ob- 
tained VCT curve has shown an upturn near T c (Ref. 



17j) in contrast to the present result. As a possible ori- 



gin of this discrepancy, let us discuss our truncation of 
the summation for Q in Eq. ([B]) which is carried out in 
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FIG. 1: (Color online) Calculated pressure to temperature 
phase diagram consisting of the corresponding VCT curve 
Tvct(P) (solid curve) in the bulk B phase or B-like phase, 
the superfluid transition curve T C (P) (dotted one), and the 
corresponding A-B transition curve Tab(P) (dashed one) of 
liquid 3 He in bulk [thin (black) curves] and an uniaxially de- 
formed aerogel [thick (red) curves]. For comparison, we also 
show the VCT curve [dashed-dotted (blue) curve] in a globally 
isotropic aerogel where the resulting T C (P) cannot be distin- 
guished from that [thick (red) dotted curve] in the anisotropic 
ones and thus, is not drawn in the figures. The parame- 
ters 8 U = 0.01 and S u — —0.003 are used for (a) compressed 
and (b) stretched aerogels, respectively, and (2-kt)^ 1 = 0.135 
(mK) was assumed in both of (a) and (b). The A- phase-core 
vortex is stable at higher temperatures than the VCT curve. 
The VCT curve in the isotropic aerogel shifts to higher tem- 
peratures than that in bulk, reflecting the reduction in the 
SC correction. As shown in (a), the A-phase-core vortex is 
stabilized by the small compression and survives even in the 
low-pressure limit, while in (b), it tends to be lost due only 
to a small stretch. 



obtaining Eq. (see the final paragraph in Sec. II). In 
the |Q| = 2 nonaxisymmetric vortex where the twofold 
symmetry is retained, we may need to include compo- 
nents with larger even \Q\ values [14J . Since such higher 
\Q\ components do not appear at f = and hardly affect 
the order parameter at the vortex core, we expect that 
the energy gain of the \Q\ = 2 double-core vortex due to 
the higher Q components mainly stems from the spacial 
variation in the order parameter far from the core which 
is relatively insensitive to temperature. Thus, inclusion 
of the higher Q components would merely shift the VCT 
curve to a higher pressure while keeping the shape of the 
VCT curve unchanged and thus, would not lead to re- 
solving the discrepancy on the VCT curve near T C (P). 
Since the problem concerning the discrepancy is not easy 



to resolve and needs further studies, we will not discuss 
the problem here. Below, we will show how the stability 
region of the A-phase-core vortex in the phase diagram 
is altered by aerogels. 

In a globally isotropic aerogel (5 U — 0) where the re- 
gion of the A-like phase is invisible in the parameter range 
used for our numerical analysis, the corresponding VCT 
curve denoted by a dashed-dotted (blue) one in Fig. 1 
shifts to higher temperatures than the bulk VCT curve, 
as a result of the reduction in the SC correction. In glob- 
ally anisotropic aerogel, the shape of the VCT curve in 
the B-like phase is drastically altered by only a small 
uniaxial deformation, depending on whether it is a com- 
pression or a stretch [see Figs. 1(a) and 1(b)]. A small 
uniaxial compression, measured by 6 U = 0.01, expands 
not only the A-like phase region but also the region of 
A-phase-core vortex while quite a small uniaxial stretch, 
measured by S u = —0.003, suppresses the region of the 
A-phase-core vortex although the A-like phase itself ex- 
pands due to the stretch. Such dependence of the A- 
phase-core vortex on the type of deformation of aerogel 
is qualitatively reasonable because, just as in the vortex- 
free A-like phase induced by a compression, 1 is aligned 
along the vortex axis in the A-phase core of the bulk B 
phase [13, LL4|. The A-phase-core vortex region in the 
stretched case becomes narrower in temperature, since 
the effect of suppression of the A-phase-core state due 
to the uniaxial stretch is so strong to overwhelm the SC 
correction just below the A-B transition curve Tab(P)- 
With increasing \5 U \, the VCT curve is extended down to 
lower temperatures in the compressed case [18| while, in 
the stretched case, the region of the A-phase-core vortex 
is suppressed and vanishes for larger \S U \ than the value 
S u = -0.006. 



IV. VORTEX CORE TRANSITION IN A-LIKE 
PHASE IN STRETCHED AEROGEL 

Next, we investigate vortex structures in A-like phases 
in aerogels deformed along the vortex axis by paying our 
attention to the deformation-induced instability of a core- 
less Mermin-Ho vortex (MHV) (Ref. [11]) in the bulk 
A phase. In uniaxially compressed aerogels, the 1 vec- 
tor tends to become parallel to the vortex axis and the 
boundary condition of MHV far from the vortex center is 
incompatible with the orientation of 1. Once the bound- 
ary condition compatible with 1 || z is chosen, the uniaxial 
compression does not affect vortex structures any longer 
and, instead, the pure phase vortex will appear. We fo- 
cus hereafter on the uniaxially stretched case in which 1 
is perpendicular to the vortex axis far from the vortex 
center. As the boundary condition, we take the follow- 
ing condition compatible with that of the dipole-locked 
MHV 
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FIG. 2: (Color online) Radial dependences of C^(r) of (a) 
an axisymmetric vortex and (b) a nonaxisymmetric one with 
\Q\ = 2 components at T=2.040(mK) and P=28 (bar) in the 
A-like phase in the uniaxially stretched aerogel with 8 U — 



<<>) 



0.03. In (a), C_^ and Cq + remain finite at 



: which 
correspond to the /3-phase component and that of MHV core, 
respectively. In (b), the components CgJ(r) are identically 
zero and other negligibly small components are not shown. 
Further, the dashed curves imply Q / components. At 
f = 0, just the two components C_j_q and C_^ remain finite, 
implying a polar core A V: i(r — 0) = a Xl z5 v ,x8i,z- 



FIG. 3: (Color online) Spacial variations in the order param- 
eter A Vt i of (a) the axisymmetric vortex and (b) the nonax- 
isymmetric one on the x axis (black curves) and y axis (red 
dashed ones). In (a), the vortex center is occupied with the 
f3 phase a x . z (x + iy) v S ltZ and the component of the dipole- 
locked MHV center a z , y S v ,z (x + iy)i- In (b), a v ,i's start 
varying around f — 10 on the y axis, corresponding to the d 
texture (see the text below), and a polar state a x ,z S v ,xSi,z is 
realized at the vortex center. The three thin solid curves in 
(a) with extremely small magnitudes denote a y , y (top), a y . x 
(middle), and a x>x (bottom). Negligibly small contributions 
on the x axis are not shown in (b). 



m = z, n = 



(9) 



with the normalization O 2 + O xv = 1. In the represen 
tation in Eq. ©, we have C_ 
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_,_ — O xy /2 with C + q = O'Iq = U z /\/2 or, equiv- 
alently, the radial 1 and d vectors. In the A-like phase 
in a uniaxially stretched aerogel, we obtain three solu- 
tions, an axisymmetric vortex, a nonaxisymmetric vortex 
with |Q| = 1 components, and a nonaxisymmetric vortex 
with \Q\ = 2 components. The \Q\ — 1 nonaxisymmet- 
ric vortex always has a higher value of the free energy 
J dfffif) than two other solutions, and thus, it will 
not be discussed hereafter. Below, we will show results 
of our calculation in the case with S u = —0.03. The re- 
sulting phase diagram is shown in Fig. [5] where a dashed 
curve denotes Tab (P) and the polar pairing state occurs 
in the narrow region sandwiched between T C (P) (dotted) 
curve and a dashed-dotted one. A solid (red) curve and 
the inset in Fig. [5] will be explained later. 

Figure [2] shows the radial dependences of Cffi (?) for 
(a) an axisymmetric vortex and (b) a nonaxisymmetric 
one with \Q\ = 2 components at T — 2.040 (mK) and 
P = 28 (bar), and Figs. |3Ia) and|3Jb) show the corre- 



sponding spacial variations in the order-parameter com- 
ponents in Ari t i(r, </>) on the x axis [solid (black) curve] 
and the y axis [dashed (red) one], where 

A Vi i(r,<p) = AA[8ri,x(a>x,xe vi 'x + ie~ vl 'a x> yy + a x ,zz)i 

+S n , z {a z , x x + i a ZtV y + e^ a z , z z)i\ (10) 

on the x (<j) — 0) and y ((f) = tt/2) axes. The axisymmetric 
vortex shown in Fig. Uta) has nonvanishing C_j_q and Cq + 
components at the vortex center r = which correspond 
to the nonunitary /3-phase component a XjZ (x + iy)r\ $z,i 
and the A-phase one a ZyX S n . z (x + iy)i in Fig. EJa), re- 
spectively. As one can see in Fig. G2a), with increasing 
?, A v ,i changes from 8 v>z [(a z , x x + a z , z z) + ia z ^ y y] t to 
8i h x(ax, z z + i a x , v y )i, indicating the continuous rotation 
of the 1 vector and the d vector from z to x. This rota- 
tional behavior in 1 and d is compatible with the texture 
in the dipole-locked MHV so that we call the axisymmet- 
ric vortex the nonunitary Mermin-Ho vortex. As one can 
infer from the sign of a z in Eqs. ([3]) and ([7]), the uniax- 
ial stretch (S u < 0) enhances a v . z corresponding to Cl . 
As temperature is lowered so that the uniaxial stretch is 



nonunitary Mermin-Ho vortex 




polar-core vortex 



FIG. 4: (Color online) Comparison of radial dependence of 
the free energy density between the axisymmetric nonunitary 
Mermin-Ho vortex and the nonaxisymmetric polar-core vor- 
tex. For both the two curves, the free energy density was 
normalized by 2nK\A\ 2 at T=2.040(mK) and P=28 (bar). 



less effective, C_j_(j(f 



0) decreases, while Cg + (f 



0) 



grows up so that A v ^ approaches that of MHV. It is 
noted that this axisymmetric vortex has a nonvanishing 
magnetic moment parallel to the vortex axis arising from 
the nonunitary /3-phase component C + q (f = 0) which 
ramains finite even at the lowest temperature in the A- 
like phase, Tab (•£")• Figure [5] (b) shows the structure of 
the nonaxisymmetric vortex, where a |Q|=2 nonaxisym- 

( — 2) 

metric part, C_ , denoted by a dashed curve remains 

nonvanishing as well as C_^ at f = 0, and the resulting 
core state becomes the nonaxisymmetric polar pairing 



one, A Vt i — a 



r, z S v , x Si,z, as shown in Fig. |3] (b). 



The radial dependence of the free-energy density /(f) 
of the polar-core vortex is also shown in Fig. 0] where 
/(f) of the nonunitary MHV is shown for comparison. 
The lower free-energy density of the polar-core vortex 
close to the vortex center competes with its bump ap- 
pearing around f = 10 in /(f). Comparing the total free 
energies f drr f(f) of the two vortices with each other, 
we find that the polar-core vortex is more stable than 
the nonunitary MHV at any temperature in the A-like 
phase. The stability of the polar-core vortex against the 
nonunitary MHV, however, depends on the vortex radius 
ffe. When f;, = 50 is assumed, the axisymmetric vortex 
has lower free energy than that of the nonaxisymmetric 
vortex in contrast the case with fj = 20. For the in- 
termediate value f b = 30, temperature dependences of 
the free energies of the two vortices at P = 28 (bar) 
are shown in the inset of Fig. [SJ where the VCT close 
to 2.02(mK) is realized. The VCT curve obtained for 
a vortex with the radius f& = 30 is denoted by a solid 
(red) curve in Fig. [5] where the polar-core vortex is sta- 
ble at higher temperatures than the solid curve. Even 
in a more realistic situation with a fixed r = R, where 
R denotes the averaged spacing between the neighboring 
vortices, the VCT should be realized because, at a fixed 
pressure, R — R/£gi,(T) becomes longer on cooling due 
to the temperature dependence of £gl(T). In this way, it 
is concluded that a discontinuous VCT fronr the nonax- 
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FIG. 5: (Color online) Calculated pressure to temperature 
phase diagram consisting of VCT (solid) curve in the A-like 
phase, T C (P) (dotted) curve, Tab(P) (dashed) curve, and the 
polar- A transition (Ref. Q) (dashed-dotted) curve in super- 
fluid 3 He in the uniaxially stretched aerogel with 8 U — —0.03. 
The polar-core vortex is stable at higher temperatures than 
the solid curve and no VCT occurs in the B-like phase. In the 
Inset, the total free energy for ft — 30 is compared between 
the nonunitary MHV (dashed curve) and the polar-core vor- 
tex (solid one) as a function of the temperature (mK). The 
VCT occurs at the temperature pointed by an arrow. 



isymmetric polar core to the axisymmertric nonunitary 
core should occur in the A-like phase. 

Finally, we comment on the origin of the bump seen 
in /(f) of the polar-core vortex in Fig. 2J As one 
can see in Fig. (3Jb), on approaching the vortex cen- 
ter along the x axis [solid (black) curves], A n ^{x) ~ 
^rj,x{io-x.yy + a-x,z z)i is almost constant in f > 4, im- 
plying that 1 vector remains nearly parallel to the radial 
direction. In contrast, on approaching the core along the 
y axis [dashed (red) curves], we have the crossover from 

-<*-r],i\y) ^ Orj,y\Q>y } x X ~r 2 CLy^z ^)i ^0 O^xyl Qx,x % t &X, z Z)i 

on sweeping through the bump in the free-energy density 
around f = 10. The spacial variation in A^j discussed 
above are schematically shown in Fig. [51 where blue and 
red arrows denote 1 and d, respectively. On approaching 
the vortex center (from right to center in Fig. [5]), the d 
vector suddenly rotates to have an alignment parallel to 
the x axis while keeping 1 in the radial direction. This 
energy cost of the resulting d texture is the origin of the 
bump seen in Fig. 3J Further, close to the vortex center, 
the orbital component perpendicular to the vortex axis 
gradually vanishes keeping d in the x direction (left in 
Fig. [S]) and the polar pairing state A^^ = a x ,z&7],x$i,z oc- 
curs at the vortex center. Therefore, we can deduce that 
a possible polar-core vortex is accompanied by a nonax- 
isymmetric d texture which, in turn, removes a singu- 
larity due to the pure phase, while keeping the 1 vector 
oriented in the radial direction. This behavior of the or- 
der parameter, implying a hybrid texture of the d vector 
and the pure phase, resembles that of the half quantum 
vortex [Tj] which might be realized in the A-like phase 
in an uniaxially compressed aerogel [7J. 
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FIG. 6: (Color online) Sketch of the 1 texture (blue arrows) 
and the d texture (red ones) in the polar-core vortex far from 
the vortex center (right) and inside the bump in /(f) (center). 
Inside the bump, order-parameter field can be expressed as 
A v> i ~ 5 VlX (a x , z z + ia x ^4>)i, indicating 1 || f and d || x, where 
a x ,<t> corresponds to a XtV for <j> — and a x , x for <j> = 7r/2 in 
Eq. (jlOp . The circle symbol and black arrow in the left figure 
expressing the close vicinity of the vortex center denote a XtZ z 
and a.x,<t><t>, respectively. In the left figure, on approaching 
the vortex center, a x ^ (i.e., the length of the black arrow) 
gradually vanishes, and at the vortex center, the polar state 
A v a = 5 n , x (a x , z z)i is realized. 



SUMMARY AND DISCUSSIONS 



A-like phase in a uniaxially stretched aerogel. The for- 
mer transition should appear even in low-pressure limit 
while it should be absent in aerogels stretched along the 
vortex axis. This conclusion might be relevant to the 
fact that no VCT has been detected in a B-like phase 
in a previous rotating measurement [191 ]. By contrast, 
it may not be easy to experimentally detect the VCT in 
the stretched A-like phase because of the random orien- 
tation of 1 due to the quenched disorder effect of aero- 
gel neglected here. Nevertheless, upon cooling the sys- 
tem through T c while rotating the aerogel including 3 Hc, 
this VCT may be seen if the lattice consisting of radial 
MHVs (Ref. [H) is formed in bulk 3 He-A under the 
similar condition because the resulting vortex lattice in 
the close vicinity of T c consists of cores overlapping with 
one another so that the order-parameter configuration 
near the core dominates over the quenched disorder ef- 
fect making the order-parameter random. An additional 
feature in vortices in anisotropic aerogels is the presence 
of a nonunitary core state which may be detected by per- 
forming a magnetization measurement. 



Possible discontinuous VCTs have been studied for 
both A-like and B-like phases in superfluid 3 He in aero- 
gels with uniaxial deformation parallel to the vortex axis. 
One VCT should be realized in the B-like phase in a uni- 
axially compressed aerogel, and a different VCT between 
a nonunitary core and a polar core may occur in the 
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